No singularities at the phase transition in the Dicke model 
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The Dicke Hamiltonian describes the simplest quantum system with atoms interacting with pho- 
tons: N two level atoms inside a perfectly reflecting cavity which allows only one electromagnetic 
mode. It has also been successfully employed to describe superconducting circuits which behave as 
artificial atoms coupled to a resonator. The system exhibits a transition to a superradiant phase 
at zero temperature. When the interaction strength reaches its critical value, both the number 
of photons and of atoms in excited states in the cavity, together with their fluctuations, exhibit 
a sudden increase from zero. Employing symmetry-adapted coherent states it is shown that these 
properties scale with the number of atoms, that their reported divergences at the critical point 
represent the limit when this number goes to infinity, and that in this limit they remain divergent 
in the superradiant phase. Analytical expressions are presented for all observables of interest, for 
any number of atoms. Comparisons with exact numerical solutions strongly support the results. 
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The quantum description of a system of N non- 
interacting two-level atoms under the action of a one 
mode electromagnetic field allowed Dicke [l| to pre- 
dict the superradiance phenomenon, where all atoms 
absorb and emit light collectively and coherently p|. 
For N=l (the Jaynes-Cummings model), the system ex- 
hibits quantum features such as collapse and revivals of 
the atomic inversion Q and squeezing of the radiation 
field [J . 

The Dicke quantum phase transition has recently been 
observed employing a superfluid gas in an optical cav- 
ity [j| . While this breakthrough was achieved using time 
dependent fields dressing the system, it has been argued 
that superconducting circuit QED systems will exhibit 
a stationary superradiant phase [||. The description of 
these systems is given in terms of the same Dicke Hamil- 
tonian Q, with minor modifications Q. The theoretical 
description is obtained by mapping the original Hamil- 
tonian to a system of two bosons, making a Holstein- 
Primakoff realization of the quasi-spin operators which 
represent the collective atomic excitation, and truncat- 
ing the infinite series up to first or second order assum- 
ing that the expectation value of the number of excited 
atoms (or the number of photons) is small compared with 
the total number of atoms N [Q-Q . This system exhibits 
a transition to a superradiant phase at zero temperature. 
When the interaction strength reaches its critical value, 
a divergence is predicted in the expectation value of the 
number of photons 0, Q > the number of atoms in excited 
states [|[ and their fluctuations 0, H • 

In this contribution we show, employing symmetry 
adapted coherent states, that these properties scale with 
the number of atoms, that their reported divergences at 
the critical point represent the limit when this number 
goes to infinity, and that in this limit they remain di- 
vergent in the superradiant phase, contrary to what is 
stated in some of the references mentioned above 



The ground state of the Dicke model is obtained using 
the tensorial product of spin coherent and Weyl states 
as a trial wavefunction. The Hamiltonian expectation 
value is calculated, and its critical points found, showing 
the existence of the quantum phase transition. The su- 
perradiant phase is described employing combinations of 
coherent states which preserve the Hamiltonian Ci sym- 
metry 10]. With them, analytical expressions are ob- 
tained for all observables of interest, which are valid for 
any number of atoms. Comparisons with exact numerical 
solutions strongly support the results. 

The Dicke Hamiltonian involves the collective interac- 
tion of N two-level atoms with energy separation equal to 
hu)A with a one mode radiation field of frequency ujf- Its 
deduction from the quantization of the electromagnetic 
field makes the following assumptions: (i) The dipole ap- 
proximation or the long wavelenght limit in such a way 
that the electric field is evaluated at the center of mass of 
the atoms and is independent of the position, (ii) Only 
the two atomic levels are interacting with the electromag- 
netic field. The lowest energy state is at least metastable 
and thus the decay to other energy levels can be ne- 
glected. The Dicke model Hamiltonian for N identical 
2-level atoms immersed in a one mode electromagnetic 
field, is given by 



H D = a^a + uj a J z H = (a* 



J-+ J 4 



(i) 



where uja = uja/uif > is given in units of the frequency 
of the field, and 7 = 7/wj? is the (adimensional) coupling 
parameter. The operators a, a', denote the one-mode 
annihilation and creation photon operators, respectively, 
J z the atomic relative population operator, and J± the 
atomic transition operators. 

The Dicke Hamiltonian has a dynamical symmetry as- 
sociated with the projectors of the symmetric and anti- 
symmetric representations of the cyclic group C'2, given 



2 



by, 



Ps = r 1 



This symmetry allows the classification of the eigenstates 
in terms of the parity of the eigenvalues A = j + m + v 
of the excitation number operator 



A= \J J 2 + 1/4- 1/2 + J z + a T a, 

with j = AT/ 2, j + m the number of atoms in their excited 
state, and v the number of photons. 

To obtain analytical expressions the variational pro- 
cedure described in [ll| is employed. We use as a trial 
state the direct product of coherent states in each sub- 
space: Heisenberg-Weyl states \a) for the photon sec- 
and SU(2) or spin states |£) for the particle 
la. i.e.. | a, Q = l a )®lC)- Due to the symme- 



tor 
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sector 

try properties of the Dickc Hamiltonian, it is convenient 
to replace m by A — j — v: Explicitly, it is given by the 
expression 



«, = Kol /2 E E rgr 

A=0 !/=max[0, A-2j] 



2j 
A- i/ 



1/2 



C A -"k)«) |i, A-iz-j) , 



(2) 



where the ket \v) is an eigenstate of the photon num- 
ber operator, \j, m) is a Dicke state with j denoting the 
eigenvalue of J 2 , and m the corresponding eigenvalue of 
J z . The trial state contains N = 2j particles distributed 
in all the possible ways between the two levels and up to 
an infinite number of photons in the cavity. The normal- 
ization is 



Ncoh = (l + ICI 2 ) exp|a| 



(3) 



The expectation value of the Dicke Hamiltonian with 
respect to |a, C), when divided by the total number of 
particles, is called the energy surface of the Dicke Hamil- 
tonian. It is a function of four variables and two param- 
eters, 



E(q, p, 9, 



— (p 2 + q 2 ) - -OJa cost 
2N K ' 2 



v^7 



q sin cos , 



(4) 



where use has been made of the harmonic oscillator re- 
alization of the field, and of the stereographic projection 
for the quasi-spin sector 
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(q + ip) 



c 



tan | - ] exp(-i0) , (5) 



(q, p) correspond to the expectation values of the quadra- 
tures of the field and {9, 0) determine a point on the 
Bloch sphere. 



The minimum and degenerate critical points are ob- 
tained by means of the catastrophe formalism (l6j . The 
minimum critical points are, for to a > 4 7 2 : 9 C = , 
q c = , and p c = ; and for < lua < 4 7 2 : 9 C — 

arccos (j^j , 0c = 0, q c = - V2N 7 1 - cos C , 
p c = . The first group contains the minimum critical 
points in the normal phase, the second one describes the 
minimum critical points in the superradiant phase. The 

/wI/2. 



phase transition occurs at lua — 4 7^:, i.e. 7c 
In dimensional units, 7 C = vSaSf/2, in agreement with 
Ref. di]. 

The associated energies, expected number of photons 
{N p h) = Wc\ 2 and of excited atoms (n e ) = + Jz) = 
y(1 — cos9 c ) are 



2 { Tc 



E norm al = -2A^7c > E super = —AT 7 \[~^ I + 1 
{^ph^normal — j {^ph} super — ^ T [ ^ 



V^ef normal — U , v^e) super — n I I 

2 \ V 7 



(6) 



These mean field expressions exactly coincide with those 
presented in and reproduce with precision the nu- 
merical results even for relatively small number of par- 
ticles. The mean field results fail badly, however, in the 
description of the fluctuations of these observables, a lim- 
itation also found in the Tavis- Cummings model, where 
the counter- rotating terms are absent (17| . The reason 
is that in the superradiant phase the state \a, breaks 
the C2 symmetry of the Hamiltonian, mixing all A val- 
ues, and having an expectation value of (q) = q c different 
from zero, implying that (Aq) 2 = 1/2. 

In order to go beyond this mean field treatment, to 
obtain a symmetry-adapted description of the ground 
state of the Dicke model in the superradiant phase, which 
would preserve the Ci symmetry of the Hamiltonian, we 
perform a projection over the A's in Eq. ((2|), which is 
equivalent to requiring that the sum be carried over A 
even or odd. This represents a generalization of the 
restoration of the Hamiltonian symmetr y in the LMG 
[l|| and Tavis- Cummings Hamiltonians {17'] - The two 
resulting orthogonal states with A even or odd \a c , Cc, ±) 



are even or odd under 



1JV-1 



A/-± = 



1 



Their normalizations are: 



(l±F(iV, 7 , 7c )). (7) 



/v 



where we have defined 
W%7c) = ^) expj-2iV7 2 h- 



(8) 
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From them analytical expressions can be found for all 
the observables of interest [10|. Those for the even and 
odd energy surfaces are 



(Aq) 2 /N 



(H) 



TV 



-3 



1±F 




(9) 



It should be noted that F(TV, 7,7 C ) — > when TV — > oo 
very rapidly. For example, F(100, 1.037 c , 7c) < 10~ 5 . 
It implies that the two states with different parity Ci 
become degenerate. Notice that the origin of this de- 
generacy is the restoration of the C2 symmetry, not its 
breaking. 

The expectation value of the population of atoms in 
the two level system is 



N 



1 - 



1 - 



1±F 



(10) 



and the expectation value of the number of photons is 
given by 



(Nphh 
TV 



1_ J±. 




1±F 



(11) 



When F —> both expressions coincide with the mean 
field ones, given in Eqs. (JBJ, as expected. 

The quadrature components of the electromagnetic 
field change the value of A in one unit, so their expecta- 
tion values with respect to the projected states are equal 
to zero. For the same reason the expectation values of 
matter observables (J x ) and( J y ) are null. The symmetry 
projection makes all the difference with respect to the 
mean field estimations, where these expectation values 
grow with the number of particles. 

The expectation values of squared quadratures of the 
one-mode electromagnetic field are given by 



TV 

(P 2 )± 
TV 



— + 2 7 2 ( 1 
2 TV ' \ 



— - 2 7 2 1 - 
2TV ' \ 









1±F 




F 




F±l 



(12) 



The expectation value of the operators J 2 and J 2 are 

(TV- 



{Jl)± 


1 


TV 2 


4TV 
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TV 2 


4TV 




In Fig. [T] the squared fluctuations of the quadrature 
(Aq) 2 at resonance (oja = 1) are presented, divided by 
the total number of particles N, as a function of the 
coupling strength 7. The thin red lines were calculated 




FIG. 1: Squared fluctuations of the first photon quadrature 
(Aq) 2 /N at resonance, oja = 1, calculated through numerical 
diagonalization of the matrix Hamiltonian for iV = 10, 20, 30 
and 50 (red thin lines ) and from Eq. 12 (black thick dashed 
line). As N increases the exact quantum dispersion ap- 
proaches the calculated through adapted states: the upper 
thin curve corresponds to TV = 10, while the lower thin line 
corresponds to N — 50. All the quantities plotted are dimen- 
sionless. (Color online) 
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FIG. 2: Squared fluctuations of the atomic transition oper- 
ator (AJ X ) 2 /N 2 at resonance, oja = 1, calculated through 
numerical diagonalization of the matrix Hamiltonian for TV = 
10, 20, 30 and 50 (red thin lines) and from Eq. 13 (black 
thick dashed line). As TV increases the exact quantum disper- 
sion approaches the calculated through adapted states: the 
upper thin curve corresponds to TV = 10, while the lower thin 
line corresponds to TV = 50. All the quantities plotted are 
dimensionless. (Color online) 



by numerical diagonalization of the Hamiltonian for N= 
20, 40, 60 and 100, while the thick green line was obtained 
employing Eq. (TT2"j) . It is clear that the fluctuations scale 
with the number of particles. 

In Fig. [5] the squared fluctuations of the atomic ob- 
servable (AJ X ) 2 at resonance (lua — 1) are presented, 
divided by the square of the total number of particles TV, 
as a function of the coupling strength 7. The thin red 
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lines were calculated by numerical diagonalization of the 
Hamiltonian for N= 20, 40, 60 and 100, while the thick 
green line was obtained employing Eq. (ITUl) . In this case 
the fluctuations scale with the square of the number of 
particles. 

It is worth to study the superradiant phase in the 
limit 7c / 7 — > 0, which is exactly solvable [19|, [HJ. The 
expectation values discussed above become: (N p h) —> 
iV 7 2 , (n e ) -»■ N/2, (Aq) 2 -> 2iV 7 2 + ±, (Ap) 2 -> 
i, £ 2 = 4 (AJ X ) 2 /AT -> JV, ^ = 4(AJ,) 2 /A -> 1. 

In this limit both the average number of photons and 
the fluctuations in the q quadrature scale with the num- 
ber of particles N and with 7 2 . The average number of 
excited atoms goes to N/2, and with it (n e )/N —> 1/2. 
It clarifies why any Taylor expansion of this quantity 
around zero, truncated to second or third order, cannot 
reproduce the asymptotic behavior. 

The squeezing in spin observables is given by £j = 
•\/4(AJi) 2 /N , where Ji are the perpendicular compo- 
nents to the average angular momentum vector (J) (2lj . 
There is squeezing when £j < 1. For the Dicke model 
(J) = (J z ) k, and the perpendicular components of the 
angular momentum are J x and J y . The squared squeez- 
ing observables £ 2 are, in the limit (n e )/N << 1, very 
close to the fluctuations of the quadratures of the atomic 
degrees of freedom studied in Q. In a similar fashion as 
the number of excited atoms, the squeezing parameter £ 2 
scales with the number of atoms. 

It is worth mentioning that Eqs. (|12|13j) show that 
(An) 2 and £ y have some squeezing around 7c , as found 
in 0. This is an effect of order 1/N, which becomes 
negligibly small when the number of atoms grows. The 
excitation energy of the first excited state is also of this 
order. These aspects will be analyzed in detail in a forth- 
coming publication. 

What is the origin of the divergences in the number of 
photons and its fluctuation at the critical value of 7 re- 
ported in 0-11? It is found in the normal phase when the 
limit N —> oo is assumed, obtaining a Hamiltonian bilin- 
ear in two kinds of bosons without memory of the number 
of particles, except for the ground state energy which has 
a constant displacement —Nuja/2. A Bogoliubov trans- 
formation allows the exact diagonalization of this Hamil- 
tonian, and the expectation values of the number of pho- 



tons, the number of atoms in excited states, and their 
fluctuations, all diverge at 7c . It agrees with the exact 
results reported above when the number of atoms N goes 
to infinity. But, in contrast to what is found in 0-0, we 
have shown that at the superradiant phase they will re- 
main divergent with N for any 7 > 7c . This finding relies 
on the use of symmetry projected SU(2)xU(l) coher- 
ent states. The expectation values obtained employing a 
truncated Hamiltonian completely disagree with the nu- 
merically exact ones at the superradiant phase because 
the assumption (n e ) « N is not valid for 7 >> 7c , and 
the Ci symmetry of the Hamiltonian is not preserved. 

In summary, we have presented analytic expressions for 
the ground and first excited states of the Dicke Hamilto- 
nian in the superradiant phase which have the appropri- 
ate C2 symmetry. From them, analytical expressions for 
all expectation values of interest were obtained, valid for 
any number of atoms. The procedure was the following: 
the expectation value of the Hamiltonian with respect to 
the tensorial product of Weyl and SU(2) coherent states 
defines an energy surface depending on phase space vari- 
ables and parameters, whose critical points were deter- 
mined. The minimum critical points yield the expres- 
sions (JSJ) which give the minimum energy of the sys- 
tem together with information about the quantum phase 
transition present in the Dicke Hamiltonian. 

At the phase transition the expectation values of the 
number of photons, the number of atoms in excited states 
and their fluctuations exhibit a sudden increase from 
zero. It was shown that these properties scale with the 
number of atoms, that their reported divergences at the 
critical point represent the limit when this number goes 
to infinity, and that in this limit they remain divergent 
in the superradiant phase. This implies that no singular 
behavior in the number of photons, the number of excited 
atoms, their fluctuations, squeezing and quadratures will 
be found around the phase transition in future experi- 
ments. The presence of entanglement [13] and chaos [l9[ 
at the phase transition in the Dicke model should also be 
analyzed from this perspective. 
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